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The method of iterated conformal maps for the study 
of Diffusion Limited Aggregates (DLA) is generalized to the 
study of Laplacian Growth Patterns and related processes. 
We emphasize the fundamental difference between these pro- 
cesses: DLA is grown serially with constant size particles, 
while Laplacian patterns are grown by advancing each bound- 
ary point in parallel, proportionally to the gradient of the 
Laplacian field. We introduce a 2-parameter family of growth 
patterns that interpolates between DLA and a discrete ver- 
sion of Laplacian growth. The ultraviolet putative finite-time 
singularities are regularized here by a minimal tip size, equiva- 
lently for all the models in this family. With this we stress that 
the difference between DLA and Laplacian growth is NOT in 
the manner of ultraviolet regularization, but rather in their 
deeply different growth rules. The fractal dimensions of the 
asymptotic patterns depend continuously on the two param- 
eters of the family, giving rise to a "phase diagram" in which 
DLA and discretized Laplacian growth are at the extreme 
ends. In particular we show that the fractal dimension of 
Laplacian growth patterns is much higher than the fractal di- 
mension of DLA, with the possibility of dimension 2 for the 
former not excluded. 



I. INTRODUCTION 

This paper had been motivated by an apparent con- 
sensus on DLA and Laplacian growth Patterns being in 
the same universality class in terms of their asymptotic 
fractal dimensions jjj . We present here a theory of these 
processes in two dimensions which clarifies the differences 
between these processes, showing in particular that their 
asymptotic fractal dimensions differ. 

Laplacian Growth Patterns are obtained when the 
boundary T of a 2-dimensional domain is grown at a 
rate proportional to the gradient of a Laplacian field P. 
Outside the domain V 2 P = 0, and each point of T is 
advanced at a rate proportional to VP In Dif- 

fusion Limited Aggregation (DLA) Q a 2-dimensional 
cluster is grown by releasing fixed size random walkers 
from infinity, allowing them to walk around until they 
hit any particle belonging to the cluster. Since the par- 
ticles are released one by one and may take arbitrarily 
long time to hit the cluster, the probability field is sta- 
tionary and in the complement of the cluster we have 
again V 2 P = 0. The boundary condition at infinity is the 
same for the two problems; in radial geometry as r — ► oo 
the flux is VP = const x f/r. Since the probability for 
a random walker to hit the boundary is again propor- 



tional to VP, one could think that in the asymptotic 
limit when the size of the particle is much smaller than 
the radius of the cluster, repeated growth events lead to 
a growth process which is similar to Laplacian Growth. 
Of course, the ultraviolet regularizations in the two pro- 
cesses were taken different; in studying Laplacian Growth 
one usually solves the problem with the boundary con- 
dition P = an where a is the surface tension and k the 
local curvature of T || . Without this (or some other) ul- 
traviolet regularization Laplacian Growth reaches a sin- 
gularity (cusps) in finite time || . In DLA the ultraviolet 
regularization is provided by the finite size of the random 
walkers. However, many researchers believed jjj that this 
difference, which for very large clusters controls only the 
smallest scales of the fractal patterns, were not relevant, 
expecting the two models to lead to the clusters with the 
same asymptotic dimensions. While we argue below that 
the difference in ultraviolet regularization is indeed not 
crucial, we maintain that the two problems are in two 
different universality classes. 

In this paper we construct a family of growth processes 
that includes DLA and a discrete version of Laplacian 
Growth as extreme members, using the same ultraviolet 
regularization. We thus expose the essential difference 
between DLA and Laplacian Growth. DLA is grown se- 
rially, with the field being updated after each particle 
growth. On the other hand all boundary points of a 
Laplacian pattern are advanced in parallel at once (pro- 
portional to VP). We show that this difference is fun- 
damental to the asymptotic dimension, putting the two 
problems in different universality classes. An announce- 
ment of these results was presented in || . 

To reach these conclusions we formulate a theory of 
Laplacian Growth patterns in terms of iterated conformal 
maps. Such a theory was successfully advanced recently 
for DLA [[7 - 10 1 , providing for an unprecedented analytic 
control of the properties of DLA |[l],[l2) ■ By generalizing 
it to Laplacian Growth patterns we can enjoy similar 
advantages, allowing us to address delicate points that 
are beyond the scope of direct numerical simulations and 
previous analytic attempts. 

In Sect. 2 we extend the iterated conformal maps ap- 
proach to parallel processes of layer-by-layer growth with 
varying local growth rates. In Sect. 3 we construct a 
2-paramctcr family of parallel growth processes that in- 
cludes DLA and the discrete Laplacian Growth as special 
(and distinct) cases. We demonstrate the relevance of 
the two parameters in determining the asymptotic frac- 
tal properties of the resulting patterns. In Sect. 4 we 
consider our algorithm for Laplacian growth and com- 
pare it to the exact dynamics without surface tension. 
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We study the correspondence between these models for 
the early dynamics (before the appearance of finite time 
singularities in the latter). In Sect. 5 we offer concluding 
remarks. 



II. ITERATED CONFORMAL MAPS FOR 
PARALLEL GROWTH PROCESSES 

The method of iterated conformal maps for DLA was 
introduced in 0. Here we present a generalization to 
parallel growth processes. We are interested in 3>( n \w) 
which conformally maps the exterior of the unit circle e 
in the mathematical w-plane onto the complement of the 
(simply-connected) cluster of n particles in the physical 
z-plane. The unit circle is mapped onto the boundary 
of the cluster. In what follows we use the fact that the 
gradient of the Laplacian field VP(z(s)) is 



1 



z(s) = & n \e w ) 



(1) 



Here s is an arc-length parametrization of the boundary. 
The map & n >(w) is constructed recursively. Suppose 
that we have already <&( n )(ui) which maps to the exterior 
of a cluster of n particles in the physical plane and we 
want to find the map <!>(™ +p ) (w) after p additional par- 
ticles were added to its boundary at once, each propor- 
tional in size to the local value of |VP|"/ 2 . To grow one 
such particle we employ the elementary map <p\fi which 
transforms the unit circle to a circle with a "bump" of 
linear size VA around the point w = e l8 . In this paper 
we employ the elementary map JtJ 



2w 



1 + w + w [ 1 + 



1 2 1 - A 



w 2 w 1 + A 



1/2' 




If we update the field after the addition of this single 
particle, then 



$(»+!) ( W ) = $(™)(^ 



(4) 



where $( n )(e' e ™ +1 ) is the point on which the (n + l)-th 
particle is grown and \/An is the size of the grown particle 

divided by the Jacobian of the map $( n ) (e l6, " +1 ) at that 
point. 

The map ^^(w) adds on a new semi-circular bump 
to the image of the unit circle under (w). The 

bumps in the z-plane simulate the accreted particles 
in the physical space formulation of the growth pro- 
cess. For the height of the bump to be proportional to 
| VP(z(s))|"/ 2 we need to choose its area proportional to 
|$(n-i)'( e ie»)|-a ( S eeEq. (0)), or 



A„ = 



A 



|$(n-l)'( e ie„)|a+2 



(5) 



With a = these rules produce a DLA cluster. Next, to 
grow p (non-overlapping) particles in parallel, we accrete 
them without updating the conformal map. In other 
words, to add a new layer of p particles when the cluster 
contains m particles, we need to choose p angles on the 
unit circle {9 m +k}k=i- At these angles we grow bumps 
which in the physical space are proportional in size to 
the gradient of the field around the m-particle cluster: 



A 



An 



ra+h 



|(j,(m)'(gi0 m+)j )|a+2 



1,2. 



(6) 



After the p particles were added, the conformal map and 
the field should be updated. In updating, we will use p 
compositions of the elementary map <ft\ t e(w). Of course, 
every composition effects a reparametrization of the unit 
circle, which has to be taken into account. To do this, 
we define a series {9 m +k}k=i according to 



Q{m)r e i6 m +k\ = «j>(m+fc-l) r e i8 m +k\ 



Next we define the conformal map used in the next layer 
growth according to 

$c™+p)(a,) = $M o 4> 9m+1 ,x m+1 0...0 to ra+ „ Ara+ » . 

(8) 

In this way we achieve the growth at the images under 
<j>(™ 1 ) f the points {0 m +k} p k=1 - To compute the 9 series 
from a given 9 series we use Eq.(||) to rewrite Eq.(0) in 
the form 



t^ m+1 ^ m+k ) ( 9 ) 



The inverse map cf> g , is given by (f>a Aui) — 
e iB (t>-\{e- l( >uj) with 



Xco 2 ± VA 2 w 4 - lu 2 [1 - (1 + AW - (1 + A)] 



1-(1 + A)w 2 



(10) 



where the positive root is taken for Re lo > and the 
negative root for Re to < 0. We stress that if we had taken 
9 n = 9 n , neglecting the effects of reparametrization, we 
would find abnormally small bumps on the tips. The 
reparametrization tends to move arcs that have to be 
mapped to fjords to regions that are mapped to tips. 
Then bumps that were supposed to grow in fjords where 
their size were normal would be pushed to tips where 
their size becomes extremely small; small bumps would 
appear where they do not belong. 

In Fig. 1 we present a schematic diagram of the parallel 
growth described above. 
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w - plane z - plane 




FIG. 1. A schematic diagram of parallel growth with four 
particles. The points 1-4 in the mathematical (u) plane stand 
for §i . . . #4, and are mapped under $' m ' to the appropriate 
images in the physical (z) plane. The sizes of the bumps were 
chosen to simulate a > 0. 

It is important to notice that on the face of it the con- 
formal map (||) appears very similar to the one obtained 
in DLA, But this is deceptive. Here we summarize 

the three major differences between DLA and the new 
growth models: 

• The distribution {0}" =1 (which is chosen uniform 
in DLA and related growth models, H) is trans- 
formed to a nonuniform (and maybe even singular) 
distribution of the angles {9}™ =1 . 

• The field by which we calculate A n (i.e. the deriva- 
tive of the conformal map) is not updated after each 
step, but only after growing a number of particles, 
and in the limit a whole layer. 

• the parameter a in Eq. can be taken to be 
different than (which is the value used to grow 
DLA clusters). In particular with a = 2 the size 
of the bump is proportional to the gradient of the 
local field, as is appropriate for Laplacian Growth. 

Note that our algorithm is not purely parallel, as the 
composition in Eq.(||) indicates. The parallel aspect is 
in using the same field to compute the values of A n in 
Eq. (g) and choosing a uniform distribution of {0}™ =1 
instead of {9}^ =1 . Note that here and below when we 
say a "uniform distribution" we mean that the series was 
created without any preference for any region of the unit 
circle. It very well may be that after avoiding overlaps 
the resulting distribution may be unevenly represented 



over the unit circle. Anyway we make use of an ordered 
series of compositions of the basic map </> to construct 
one layer. In the next section we will show that the serial 
aspect of the layer growth is not important in terms of 
the asymptotic fractal dimension of the clusters. In other 
words, the order of placing the bumps is not relevant as 
long as the same field is used as in Eq. (||). 

The details of the algorithm, including how to choose 
the series {O m +kYk=i to av °id overlaps are presented in 
Appendix A. 

III. TWO-PARAMETER FAMILY OF GROWTH 
PROCESSES 

Evidently, a discretized Laplacian Growth calls for 
choosing the series {O m +kYk=i sucn as to nave full cover- 
age of the unit circle (implying the same for the boundary 
r). On the other hand DLA calls for growing a single par- 
ticle before updating the field. Since it was shown that 
in DLA growth A„ decreases on the average when n in- 
creases, in the limit of large clusters DLA is consistent 
with vanishingly small coverage of the unit circle. To in- 
terpolate between these two cases we introduce a param- 
eter that serves to distinguish one growth model from the 
other, giving us a 1-parameter control (the other parame- 
ter is a). This parameter is the degree of coverage. Since 
the area covered by the pre-image of the n-th particle on 
the unit circle is approximately we introduce the 

parameter 

I I 

C = ~Y V^n+k ■ (11) 

(In the Appendix A we show how to measure the coverage 
exactly) . Since this is the fraction of the unit circle which 
is covered in each layer, the limit of Laplacian Growth is 
obtained with C = 1. DLA is asymptotically consistent 
with C — 0. Of course, the two models differ also in 
the size of the growing bumps, with DLA having fixed 
size particles, (a = in Eq.(||), and Laplacian Growth 
having particles proportional to VP (a = 2 in Eq.(||)). 
Together with C we have a two parameter control on 
the parallel growth dynamics, with DLA and Laplacian 
Growth occupying two corners of the a,C plane, at the 
points (0,0) and (2,1) respectively. 

Needless to say, with our partially serial growth within 
the layer, we introduced an extra freedom which is the 
order of placement of the bumps on the unit circle. In 
order for the model to have a physical meaning (i.e. to 
simulate true Laplacian growth), it must not depend on 
the specific itinerary of {fti}^m+i used to cover the unit 
circle as long as it is uniformly distributed on the unit 
circle. We will show first that this extra freedom has no 
consequence with regards to the asymptotic dimension of 
the resulting cluster. For fixed a the dimension depends 
only on the value of C. To demonstrate this, we will 
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consider various itineraries to achieve a uniform coverage 
C. 

One way is to construct the "golden mean trajectory" 
dm+k+l — Om+k + 27rp 

where p = (\fb — l)/2. At each step we check whether 
the newly grown bump may overlap a previous one in the 
layer. If it does, this growth step is skipped and the orbit 
continues until a fraction C is covered. The first bump of 
the next layer is grown at a random position in order to 
eliminate correlations induced by the arbitrary itinerary 
chosen to grow the previous layer. Another method is 
random choices of 9 m +k with the same rule of skipping 
overlaps. A third method is what was termed in |H] the 
"period doubling" itinerary 

% = 0, 9 2 n +k = k + ^ , < k < 2" , n > 

(12) 

In all these methods we cannot reach C = 1, since there 
are gaps left between the bumps. These are excluded 
from further growth in the present layer since their sizes 
are smaller than the corresponding value of yAn as cal- 
culated in the middle of the gap. We can estimate the 
maximal value of C to be of the order of 0.65. Never- 
theless, to be an acceptable model of parallel growth the 
fractal dimension of the resulting cluster should be in- 
variant to the itinerary. This invariance is demonstrated 
below. In the next section we treat the case C = 1 by 
constructing an ordered series with extra care. In the 
rest of this section we demonstrate the irrelevance of the 
itinerary. 

With the present technique it is straightforward to de- 
termine the dimension of the resulting cluster. The con- 
formal map <$>( n \uj) admits a Laurent expansion 

$W (w) = F[ n) u + F (n) + ^+ .... (13) 

The coefficient of the linear term is the Laplace radius, 
and was shown to scale like 

p (n) _ S l /D _ (w) 

where S is the area of the cluster, 

n 

S = £> |$0-D'(e*')| a ■ (15) 

Note that for a = this and equation (||) imply that S — 
jiXq. Indeed for a = this estimate had been carefully 
analyzed and substantiated (up to a factor) in Q. On 

(n) 

the other hand F± is given analytically by 

n 

p[ n] = n votm~) , (i6) 

k=l 



and therefore can be determined very accurately. 

In order to achieve comparable growth rates for differ- 
ent layers we inflated Ao in Eq. @ according to Ao — ► m\ 
in the layer composed of p particles {m + fc}? =1 . The ex- 
act form of inflation is not important; we introduce it 
simply to oppose the slowing down due to the decrease 
of (A n ) with n. In Fig. 2 we show F% of clusters grown 
by choosing the three different itineraries discussed above 
to produce the layers and for two values of C. We note 
that the curves superpose for the 3 different clusters with 
the same value of C; for different values of C a different 
behavior of F\ is manifested. 

1000 e a 




10 100 1000 10000 
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FIG. 2. Log-log plots of Fi vs. S of six individual clusters 
with a = 2, using 3 different itineraries for layer construction, 
with two values of C. C = 0.3 (upper group) and C = 0.5 
(lower group). Here we use the golden-mean, random and the 
period doubling itineraries (see Ref. [9]). 

We conclude that the dimension (determined by the 
asymptotic behavior of F\ vs. S) does not depend on the 
itinerary used to form the layers but on C only. To un- 
derstand this further we note that changing C for a fixed 
value of a is equivalent to changing the growth probabil- 
ity. For example, if we perform growth with C — 1 and 
a = we generate compact clusters, since we grow fixed 
size particles with a uniform measure. C = and a = is 
DLA even if we choose a fixed, or even a growing number 
of particles in each layer, as long as C — * asymptoti- 
cally. In Fig. 3 we demonstrate this fact by simulating 
growth with 1 particle in each layer (DLA), and growth 
with random addition of particles to each layer until the 
first overlap. In both cases C = 0, as the number of parti- 
cles in each layer grows slower than the available number 
of sites, and see Sect. 5 for a proof of this statement. 
Accordingly the dimension is invariant, despite the fact 
that the number of particles in each layer increases to 
infinity with the cluster size. 



4 




FIG. 3. (a) Regularly grown DLA (1 particle in a layer), 
(b) Cluster grown by randomly attaching bumps in each layer 
until the first overlap. Both clusters contain 10 000 particles 



In Fig. 4 we show three fractal patterns grown with 
three different values of C > 0. We draw the reader's 
attention to the fact that drawing cluster like the one in 
panel c is not entirely trivial. Simply mapping the unit 
circle will not work since many of the fjords will be lost. 
In fact, in appendix |c| we develop a reliable and effective 
method to produce the border of the fractal cluster. 



a % 




b 







FIG. 4. Patterns grown with a = 2 and 3 different values of 
C by using the golden-mean itinerary: a) C — 0.1, b) C = 0.3, 
c) C = 0.5. 



Even a cursory observation of these patterns should 
convince the reader that the dimension of these patterns 
grows upon increasing C. For a quantitative determina- 
tion of the dimension we averaged F\ of clusters produced 
by the golden mean itinerary, each with another random 
initial angle in each layer. Plots of the averages (Fi) for 
3 values of C are presented in Fig. 5. 
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FIG. 5. Linear regressions of log- log plots of (Fi) vs. S for 
a = 2 and 3 values of C: 0.1 (solid line), 0.3 (dotted) and 0.6 
(dashed). The slopes of the curves imply dimensions D=1.37, 
D=1.75 and D=1.85 respectively. The averages are taken over 
20 clusters. 

We conclude that the dimension of the growth pattern 
increases monotonically with C, with D ss 1.85 when C = 
0.6. 

In Fig. 6 we present the a, C "phase diagram" which 
results from calculations for a variety of values of C and 



into the fjords not allowing them to hit the tips only as 
is highly probable. Also decreasing a to a — increases 
the dimension, since we grow equal size particles into the 
fjords, whereas increasing a reduces the size of particles 
added to fjords and increases the size of particles that 
accrete onto tips. In particular it is obvious that DLA 
and our discretized Laplacian Growth cannot have the 
same dimensions, putting them in different universality 
classes. In particular the dimension D = 1.85 obtained 
for a = 2 and C = 0.6 is a lower bound for the dimension 
of Laplacian Growth patterns. This is because the di- 
mension increase with C and C — 1 for Laplacian Growth 
patterns. In the next section we present evidence, by 
reconsidering C = 1, that the crucial difference between 
DLA and Laplacian Growth is not in the the discretiza- 
tion or in the ultraviolet regularization, but rather stems 
from the different values of a and C. 

Before turning to models with C = 1 we note in pass- 
ing that the present family of models warrants further 
study on its own right, independently of the relation 
between DLA and Laplacian Growth. The wealth of 
growth structures seen in electro-deposition, dielectric 
breakdown models, and bacterial colony growth |l| may 
very well justify 2-parameter families of models. The 
present one is not less physical than any other that had 
been studied so far in the literature, but it enjoys the 
benefit of easily obtained conformal formulation. 
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FIG. 6. "Phase diagram" in which the fractal dimension D 
is displayed for selected values of the parameters C and a. 

The conclusion from these calculations is that the frac- 
tal dimension of the clusters depends continuously on 
the parameters, growing monotonically upon decreasing 
a or increasing C. It is quite obvious why increasing C 
should increase the dimension, we simply force particles 



IV. DISCRETE VERSUS CONTINUOUS 
LAPLACIAN GROWTH 

Continuous Laplacian growth without surface tension 
has been studied using dynamics of conformal maps in 
^J,||. The dynamical equation for the conformal map 
reads 



Re{u&(uj,t)$ t (ui,t)} = 1 



(17) 



As is well known, the solutions of this equation generate 
finite time singularities from smooth initial data. The 
simplest example is the initial condition 



$(w,0) = Fi(0)u;- 



(18) 



The number of Laurent coefficients is preserved by 
Eq.@ with 



F-2{t) 



const 



(19) 



The finite time singularity is seen from the analytic result 
(writing F = Fj(0)) 



Fi(t) 



2F-2 



[F? - 4F? 2 ]2 _ 8tF* 



4F. 



F_ 2 (t) = — - [F? - J[F? 4^ 2 ]2 _ 8tF 2 



(20) 
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At t -> t c = [Ff - 4Fl 2 ] 2 /8F* 2 , Fi(t)/F_ 2 (t) -> 2 and 
a cusp is developed at the images of 1, exp(^p), and 
exp(^), see Fig. 7 . This simple example motivated at- 
tempts to understand the role of surface tension as an ul- 
traviolet regularization, see @. We will use this result to 
study further the correspondence between our discretized 
Laplacian growth and the continuous counterpart, and 
to solidify the fundamental difference between the latter 
process and DLA. 




FIG. 7. Evolution of cusps starting from smooth initial 
conditions, with i*i(0) = 1, F_ 2 (0) = 0.24. Curves are shown 
at initial time t = 0, an intermediate time and at the critical 
time t = t c . 

As explained in Sec. 3, reaching C = 1 is impossi- 
ble with any of the itineraries discussed above. We 
can achieve this limit by growing in an ordered fashion, 
adding bumps in a controlled manner, precisely such as 
to glue one branch cut to its neighboring one. How to 
do this while imposing the appropriate symmetries is ex- 
plained in detail in Appendices A and B. We discover 
that the growth patterns constructed in this way tend to 
fractalize rapidly due to the existence of the branch cuts, 
in agreement with our statement above that the result of 
our process is a faithful lower bound to the dimension of 
continuous Laplacian Growth. An example of the pat- 
terns grown by our discretized process from the initial 
conditions (18) is shown in Fig. 8. 
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FIG. 8. Cluster grown with C — 1 starting from the same 
initial conditions as in Fig. 7. Notice that the branch cuts lead 
to spurious fractalization of the smooth envelope. 

The idea of this section is to isolate the effects of the 
parameter C and a from effects of discretization and ul- 
traviolet regularization. To this end we eliminate the 
instabilities caused by the bumpiness by keeping track of 
the two Laurent coefficients F\ and F_ 2 . We start with 
the initial conditions Fi(0)a-> + F!_2(0)/a; 2 . Every layer is 
then grown by our algorithm with a chosen values of C 
and a, computing the new values of F± and FI2, using 
the analytic formulae presented in || . Discarding all the 
other Laurent coefficients we have an updated conformal 



map in the form F^'to + F^ 2 /lu 2 . 

We find the results of this exercise quite revealing. In 

(n) (n) 

Fig. 9 we show the computed values of F± and F_ 2 and 



? (n) , 



the ratio (|19|) for C — 1 and a — 2, together with other 
values of these parameters. For C — 1 and a — 2 the so- 
lution approximates rather closely the exact results up to 
the creation of the finite time singularity, with large devi- 
ations appearing only when the tip radius of curvature is 
of the order of Ao . The degree of approximation improves 
when Ao is reduced. On the other hand, the same proce- 
dure with other values of C or a deviates from the exact 
results immediately, with the degree of deviation being 
monotonic in the difference in values of C from unity and 
of a from 2. 
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FIG. 9. Ff/F-i as a function of S for the smooth process 
described in the text, and for the following values of C and a: 
circles, C = 1 and a = 2; squares, C = 1 and a = 0; diamonds, 
C = 1 and a = 4; triangles, C = 0.5 and a — 2. The solid 
line represents the initial conditions which remain constant, 
Eq.(19). 




FIG. 10. Fi and F-2 for the smooth growth process de- 
scribed in the text with C = 1 and three different values of 
Ao. Circles Ao = 10~ 4 , triangles Ao = 5 * 10 -5 and squares 
Ao = 10 -5 . The solid line results from solving (17) with the 
same initial conditions. 

It is not uninteresting to note the similarity between 
the effect of the finite size particles and surface tension. 
This is demonstrated in Fig. 10. The deviation from the 
analytic solution depends on Aq. The smaller the latter 



is, the deeper we go into the cusp formation, and the 
closer we get to the singularity time t c . We estimate the 
time of deviation by comparing the radius of curvature 
to the physical size of our particle at the tip. This means 
that at the tip 



Ao 



1 



|$'(tip)| 2 K 2 



(21) 



where K is the curvature at the tip. The RHS vanishes 
when t — > t c , inhibited here by the value of Ao- The time 
of deviation is therefore when Ao = |$'(tip)| 2 (i)/ft 2 (i). 
We can compute the quantities involved analytically: 



$'(ti P ) = Fi — 2F_ 2 , 
= F\ + 4F_2 



(22) 
(23) 



Accordingly we can estimate the time of deviation and 
compare it with the numerics. The agreement is excel- 
lent. 

At this point it is worthwhile to reexamine the con- 
sensus formed in favor of DLA and Laplacian Growth 
being in the same universality class. Superficially one 
could say that in DLA the update of the harmonic mea- 
sure after each particle is not so crucial, since the effect 
of such an update is relatively local [ ^5| . Thus it may 
just work that a full layer of particles would be added 
to the cluster before major interaction between different 
growth events takes place. However this view is com- 
pletely wrong. An incoming random walker lands on top 
of a previously attached one very often. To see this, con- 
sider how many angels {0j} can be chosen randomly on 
the unit circle before the first overlap between bumps of 
linear sizes ej — y/ \ n (e* 9j ). To get the order of magni- 
tude take €j = e = (\/A^). The average number of times 
that we can choose randomly an angle before the first 
overlap is Af(e) ~ ^j. The Length of the unit circle that 
is covered at that time by the already chosen bumps is 
C(e) = eAf(e) ~ y/e. It was shown in § that for DLA 
(A n ) ~ j[, so that e ~ implying Af(n) ~ n 1 / 4 . No- 
tice that this result means in particular that for a DLA 
cluster of 1 million particles only less than 50 random 
walkers can be attached before two of them will arrive at 
the same site! Moreover, C(n) ~ -^jj — > for n — > oo, 
which means that as the DLA cluster grows, our cover- 
age parameter C goes to zero, rather than to unity where 
Laplacian Growth is. Taking spatial fluctuations of A„ 
into account may change the exact exponents but not 
the qualitative result. This argument clarifies the pro- 
found difference between growing a whole layer simul- 
taneously and particle-by-particle. Note however that 
DLA is NOT the C — > limit of our 1-parameter family 
because a = in Eq. (g) for DLA and a — 2 in Eq. 
(||) for Laplacian patterns. We will now show that if we 
eliminate the basic instability that stems from particles 
landing on each other then DLA and Laplacian growth 
coincide. To do so we start again with the initial con- 
ditions F\{0)u) + F_ 2 (0)/u 2 , grow one particle with the 
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DLA rules, compute the new value of F\ and -F_2, and 
use the new map uj + / uj 2 as "initial conditions" 
for an additional particle growth. The results of this pro- 
cess are shown in Fig.ll, which is now indistinguishable 
from Laplacian growth with C = 1. In fact, when the in- 
stability produced by particles landing one on top of the 
other is eliminated, the choice of the bumps according to 
the harmonic measure simulates the growth of a layer, as 
was expected by many researchers. The presence of the 
instability which is intimately linked to the DLA growth 
rules makes it fundamentally different from the paral- 
lel layer growth of Laplacian dynamics. We believe that 
with this discussion wc offered a conclusive evidence for 
the fundamental difference between DLA and Laplacian 
Growth. 



3.0 




0.0 1 ■ 1 ■ 1 ■ 1 

0.0 0.5 1.0 1.5 

S 

FIG. 11. Growth Patterns as in Fig. 10, starting from the 
same initial conditions but growing particle by particle ac- 
cording to the DLA rules, preserving only F\ and F_2. The 
inset shows a zoom of Fi close to the time singularity. 

V. CONCLUSIONS AND REMARKS 

We have introduced a 2-parameter family of growth 
patterns with the aim of clearly separating DLA from 
Laplacian Growth. We explained how to grow in parallel, 
taking care of the delicate issue of reparametrization. For 
the latter issue we needed the inverse map as explained 
in Sect. 2. The tools developed to study and control 
the reparametrization are further employed to develop 
symmetry preserving growth algorithms (Appendix B) 
and efficient mehtods to construct the interface of fractal 
clusters (Apeendix C). We argued that the parameters C 
and a are relevant for the asymptotic dynamics, whether 
the order of placing the bumps is not. The dimensions 
of the resulting growth patterns were shown to depend 
continuously on the two parameters. Besides providing 



us with a new model which is interesting by itself, we 
could reach the following main conclusions: 

• DLA and Laplacian Growth are not in the same 
universality class. 

• The dimension of Laplacian Growth patterns had 
been bounded from below by 1.85. We do not have 
a sharp estimate of this dimension, and cannot ex- 
clude D = 2. 

• The difference between DLA and Laplacian Growth 
models is not in the ultraviolet regularization. We 
explained that the deep difference is between the se- 
rial and parallel growth events, leading to increased 
tendency to form spikes in DLA. 

In future work it may be worthwhile to attempt to 
find a sharper estimate of the dimension of Laplacian 
Growth patterns. It seems also worthwhile to study the 
connection of the present model to models with noise- 
reduction, and to further understand how to interpret the 
rich phenomenology of electro-deposition and bacterial 
colony growth. 
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APPENDIX A: DETAILS OF THE ALGORITHM 

This appendix consists of three parts. In the first we 
explain how the absence of overlaps between grown par- 
ticles can be defined in terms of the conformal map. The 
second part of this appendix is dedicated to a detailed 
description of the algorithm that was introduced in Sect. 
2. In the last Subsection we explain the algorithm used 
to achieve C = 1 at each layer. 

1. overlaps in terms of iterated conformal maps 

Suppose that the first particle in a new layer is the 
(m+ l)-th particle in the growth process, and that there 
are no overlaps between the first k particles grown in this 
layer. 

In order to express this in terms of the iterated confor- 
mal map formalism, let us make the following definitions: 

• oj^' L are the two branch points of the map 4>§ \ > 
denoted as "right" and "left" respectively. 
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u>n' L are the two branch points of the map 4>e n ,\ n , 
(which we denote in the sequel as 4> n for brevity). 

Let us further denote: 



(Al) 



Note that |/3„ — |/ (27r) is the fraction of the unit circle 
covered by the particle. 

The angles j3^ ,L and aig\Cj^ L ] are connected to each 
other in a similar manner to the way n is connected to 
On <§■ 



C-l°---°^m+l(^n) • 



(A2) 



Notice that the inverse function is analytic on the 
unit circle only outside the arc 

If the k-th particle does not overlap any of the previ- 
ously grown particles of the layer, then the three points 

are all in the image of the unit circle under $( m ). In 
other words, Eqs. ([)]) and ( |A2| ) (for n = m + fc) are 
solvable. 

Since j/)" 1 is analytic on the unit circle only outside the 
arc [j3% , Pn'} > the existence of the following set of k — 1 
conditions is necessary and sufficient for the solvability 
of Eq. (|). 



2. Define 9 m+1 = m +i- 

3. Calculate \ m +i from Eq. (^|), using the derivative 
of <J>(™) at the point e l9 '"+ 1 . 

4. Calculate /J^'-fi from Eq. (Al) and store them. 

5. Let d H/^+i-^+il/^) 

For k > 1: 

6. Calculate A m+ fc by the derivative of $( m ) at the 
point e l6lm+fc , and find the appropriate branch 
points uj m+k . 

7. Check the 2(k — 1) conditions, given in Eq. ( [A3] ) 
upon replacing 9 m +k by arg[a}^'f fe ] on the LHS of 
Eq. (A3). If any of them is violated (which means 
that the fc-th particle overlaps one of the former 
k — 1 particles in the layer), choose another 9 m+k 
and repeat from stage 6. 



3. Solve Eqs. (|) and © to find 9 m+k , p* +k , /3£ +fe , 
and store them. 



9. Let C k = C fc _! + \p* 



m-j-k 



L 



m-j-k 



1/(2*0 



10. After a series of p "good" angles {9 m +k} I k= i was 
found, such that all the p(p — 1) solvability condi- 
tions resulting from p iterations of Eq. (AS) are 
fulfilled, update the conformal map according to 
Eq. (§) and C = C p . 

In our simulation p is not constant, but is determined by 
the value of C that we want to achieve in each layer. 



+fe-U 



(A3) 



Two similar sets of k — 1 conditions each are obtained 
by substituting argp^^J instead of 9 m +k on the LHS of 
([A3]), and their existence is necessary and sufficient for 
the solvability of Eq. (|A|). 

A failure of any of this conditions means that the fc-th 
particle is overlapping at least one of the previous k — 1 
particles. It is clear that if the two edge points of the 
particle are on the boundary of the m-particles cluster, 
so must be its tip (except very rare fill-up events S that 
we can safely neglect here). Therefore the existence of 
the last two sets of k — 1 conditions each is sufficient for 



3. A full coverage (C = 1) growth algorithm 

To reach C = 1, we construct recursively a series of 
consecutive angles {djYpJZ+i sucn that the left branch 
cut of the jth particle coincides with the right branch cut 
of the (j + l)th particle. This reads 



Pf+i = arg^+iH+i)] = arg[^] 



(A4) 



(A5) 



the existence of the k — 1 conditions (A3) 



2. Growth algorithms for C < 1 

The algorithm for growing one layer of p particles on 
the cluster made up of m particles (given $(™)) is defined 
as follows: 

1. Choose a series {0m+k}^ =1 uniformly distributed 
on the interval [0, 2ir}. 



Given a pair (6j, Xj) (and hence uf' R and (3j' R ) we have 
to choose 0j + i such that the value of which is deter- 
mined by 6j+i and the value of Xj+i computed at 9j+i 
coincides with the previously computed arg[w^]. Numer- 
ically this is obtained as follows. We start with 9j+i far 
enough from 9j. Then, using Eqs.(|]),(||),(^l]), we calcu- 
late the appropriate values of 0j+i, Xj+\ and Pf + \- This 
process is repeated until a value of 0j+± is found such 
that < @f +1 - arg[^ L ] < 0.01 
the whole circle is covered. 



Xj. We proceed until 
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APPENDIX B: IMPOSING SYMMETRIES ON 
THE ITERATION SCHEME 

In this appendix we explain how to use iterations of 
conformal maps to describe growth in geometries less 
symmetric than the radial. In addition we show how 
to preserve symmetries of the continuous Laplacian dy- 
namics along the iterations. The basic idea will be 
demonstrated through the important example of growth 
in channel geometry, and straightforwardly employed to 
growth from initial conditions with reflection symmetry 
or n-fold symmetry in radial geometry. 

1. Growth in a channel 

The simplest symmetry that is preserved in the it- 
erations scheme is 27r-periodicity. Clearly, 4>e,\(^) — 
<f>8,\(e l< ' 27r+ ^)- Therefore, if the initial conditions (i.e. 
$(°)) have the property 

^°\e i ^ +e ) = ^(e ie ) + L (Bl) 

where L is the channel width, then will have this 
property for any n > 0. The simplest $^°- ) that has 
the periodicity property is of course $ < - - ) (w) — — log(w) 
(L = 1) which describes a growth starting from a flat 
curve. Notice that the boundary conditions of the Lapla- 
cian held VP = ng^ij a t infinity will be automatically 

changed from VP ~ to VP ~ const x. 

Suppose now that we want to describe a growth in a 
channel with no-flow boundary conditions at the walls. 
This means that the Laplace problem has to be solved at 
each stage with the extra boundary conditions that the 
two walls y — and y — L are streamlines of the scalar 
field P (i.e. ^| y =o,L = 0). Pre-images of streamlines 
of P in the physical plane are rays (arg(w) = const) 
in the mathematical plane. Therefore, imposing no-flow 
boundary conditions at the walls amounts to demanding 
that the two rays arg(cj) = ±e (e — * 0) are mapped under 
to the walls y = and y = L respectively, for every 

n. 

Clearly, the elementary map <frg y \((j) does not have this 
property. Except for 8 = 0,ir the ray arg(w) = is 
mapped to a curved line in the z-plane. Therefore, the 
appropriate boundary conditions at the walls are not re- 
spected by the iteration process. 

We can overcome this difficulty in an analogous way 
to the image method used in electrostatics. Given initial 
conditions defined by some < I >< - * ) we construct our 
by: 

arg[c<j] < 7r 

$ (0) (u) = l> (0) (27r-w)+Z arg[w]>7r (B2) 

Under each half of the unit circle is mapped to an- 
other copy of the original interface with reflection symme- 
try around the real axis (arg[w] = 0,7r). The pre-image 



of the two walls y — and y = L under 1 are the 

rays arg(w) = + and arg(w) = n (or arg(w) = 0~ and 
arg(w) = 7r) respectively. 

Now we construct an elementary conformal function 
that maps the rays arg(w) = 0, n onto themselves. This 
can be achieved by choosing the elementary map to be 

<Pe,\ ° 0e,AM , (B3) 

such that the image of the unit circle will have the real 
line as a symmetry axis. This is shown schematically in 
Fig. 12. 




FIG. 12. Iterative conformal function that maps at each 
stage n the unit circle and the real axis in the mathemati- 
cal plane to the evolving interface and the channel walls in 
the physical plane, respectively. The two rays arg[oj] = 0, n 
are mapped under <j>e,\ ° 4>s,\ to themselves, and under the 
operation of $ ( - ' to the walls y — 0, L. 

Since the rays arg(w) = 0, tt are mapped to the walls 
under <E>(°) they will be mapped to the walls under 
defined by 

$ (,l) (w) = $ (0) O (f>g lM O (f) SiM 0...0 (j) 9niXn O 4>0 ntXn { W ) 

(B4) 

Naively, one may think that 9 = —9. However, con- 
structing the symmetric map as a composition of two 
non-symmetric maps leads to some complication. In or- 
der to have a symmetric image of the unit circle, one 
would like to have the second bump in the image of the 
unit circle to be located exactly symmetrically to the first 
bump: 

<Pe,x(e iS ) = e~ ie ■ (B5) 
Eq. ([B5| ) implies choosing 9 according to 

6 = ^14,^-^)]. (B6) 

The difference \0 — (—9)\ becomes smaller with A, and is 
zero at the points 9 — 0, it, ±7r/2 for every value of A. 
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2. Preserving symmetries of the Laplacian dynamics 

The simple technique that was developed in the pre- 
vious subsection can be generalized for cases in which a 
symmetry of the Shraiman-Bensimon Eq. Jr7| ) is known 
for specific initial conditions and we want to preserve it 
upon using iterations of conformal maps. 



a. reflection symmetry in radial geometry 

Suppose that the initial interface has a reflection sym- 
metry with respect to some axis. Without loss of gen- 
erality we can take the symmetry axis to be the x-axis, 
which is the image under <&W of the real axis in the 
mathematical plane. Then: 



(B7) 



It is easy to prove that this symmetry will be preserved 
under the Shraiman-Bensimon dynamics. 

In order to respect this symmetry in ou r iterative 
scheme we use again the elementary map (B3) which has 
reflection symmetry with respect to the real axis. Thus, 
$ (n) (w), defined by Eq. @ with that has the 
property ( |B7| ) will preserve reflection symmetry. 



b. n-fold symmetry in radial geometry 

The Shraiman-Bensimon Equations preserve also 71- 
fold symmetry. Therefore, if the initial interface, defined 
by has this symmetry, so should do $("). For sim- 
plicity let us consider 3-fold symmetry of the form: 



In order for this symmetry to be preserved, the elemen- 
tary map must be 3-fold symmetric as well. Following 
the discussion in the first part of this appendix this can 
be achieved by choosing the elementary map to be 



<t>e,\ o (/>§,\(/>§ A M > ( B9 ) 



where 



= arg[0-i(e^0)] 

^argferio^e 4 ^)]. (BIO) 

The evolution equation for &( n > now reads 

= $(o) o ^ 6lM o ^ uXi o 5iiAi o . . . o 

0...0 tj>e n ,x n o (t>g n ,x n ° ( /'e„,A„( w ) • ( Bu ) 
The extension to higher symmetries is straightforward. 



APPENDIX C: CONSTRUCTING AN OUTLINE 
FROM BRANCH POINTS 

The common method || to produce the outline of n- 
particles cluster constructed by the iterated conformal 
map technique is to sample the unit circle at K an- 
gles {Ok}k=i an d ^0 plot their images under the map 
{$( n )(e l9fc )}^L 1 . This simple method is problematic since 
a uniform series {Ok] will sample the tips much more than 
the fjords, and thus in order to have a reasonable image 
of the fjords (which are the major part of the fractal 
cluster), a huge number K n has to be used. Since 
calculation of each image point & n )(e l9k ) calls for 0(n 2 ) 
operations, this turns out to be a very inefficient method. 

Here we propose an algorithm of 0(n 2 ) complexity to 
produce an exhaustive real-space image of the whole clus- 
ter. The key idea is to focus attention on the edge points 
of the particles, which are the images of the branch points 
of the map on the unit circle. Each growing particle 
adds on two new branch points to the evolving map and 
may remove some old ones due to overlaps (see discus- 
sion in Appendix A). Therefore, the number of "exposed" 
branch points of is bounded by 2n. Let us denote 
these points {uk R ' L }k=i- An exposed branch point u>k R ' L 
was added to the conformal map by the fc-th growing par- 
ticle, and since this particle was not overlapped by any 
of the next n — k particles it remains as a branch point 
of the map Nevertheless, the reparametrization of 

the unit circle induced by the following n — k iterations 
changes the pre-image of each branch point from u>k R ' L 
to u>k,n ■ The connection between u>k R ' L and ajk,-, 
is given, similarly to Eq. (Q) by: 



R,L 



*< n >(a£f) , (CI) 



(B8) which can be simplified to 



R,L 
} k,n 



l,A„ 



The solvability of Eq. 



(C2) determines whether the ap- 
re k-th particle remains ex- 
k particles. 



propriate edge point of 
posed under the addition of the next n 
Checking the solvability conditions and calculati ng the 
reparametrized branch points u>k,n R ' L from Eq. (C2) is 
performed in the same way as in Appendix A, and it con- 
sists of 0([k — n] 2 ) operations. The total complexity of 
the algorithm is therefor 0(n 2 ). 
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